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Abstract. The paper formulates Maxwell's equations in 4-dimensional Euclidean 
space by embedding the electromagnetic vector potential in the frame vector go. 
Relativistic electrodynamics is the first problem tackled; in spite of using a geometry 
radically different from that of special relativity, the paper derives relativistic 
electrodynamics from space curvature. Maxwell's equations are then formulated and 
solved for free space providing solutions which rotate the vector potential on a plane; 
these solutions are shown equivalent to the usual spacetime formulation and arc then 
discussed in terms of the hypersphere model of the Universe recently proposed by the 
author. 
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1. Introduction 

In recent years I've been proposing an extension of optics to 4- dimensional space, 
naturally calling this discipline ^-dimensional optics (4D0). Initially this was proposed 
as an alternative formulation for relativistic problems; two papers that use this view and 
provide a good introduction are Almeida 0, 01 • The latter of these references has some 
flaws in its approach to electromagnetism which will be corrected here; it also suggested 
a metric for vacuum which was superseded by a more appropriate one in |^. In this very 
recent work I made some cosmological predictions arising from a development of 4D0, 
which had so far been impossible with a general relativity approach, showing that dark 
matter is really unnecessary for the explanation of observations if one accepts 4-space as 
being Euclidean with coordinate x° the radius of an hypersphere. The latter work left 
unexplained the fact that photons were constrained to great circles in a 4-dimensional 
hypersphere rather than following straight line geodesies. 

The explanation of photon behaviour calls for a full exposition of Maxwell's 
equations in 4D0 context, which has not yet been done in a formal way. The 
present paper is a presentation of electromagnetism in 4D0 space, introducing the 
electromagnetic vector potential as part of the space frame and deriving Maxwell's 
equations in a natural way. The solution of Maxwell's equations leading to 
electromagnetic waves is then shown to be bound to great circles on the 4D hypersphere, 
thus solving the difficulty in Almeida j^. 

In the exposition I will make full use of an extraordinary and little known 
mathematical tool called geometric (Clifford) algebra, which received an important 
thrust with the works of David Hestenes ^. A good introduction to geometric algebra 
can be found in Gull et al. and in the following paragraphs I will use the notation and 
conventions of the latter. Expressing Maxwell's equations in the formalism of geometric 
algebra is not new; the Cambridge Group responsible for the reference above uses this 
approach in one of their courses p] for the relativistic formulation of those equations, 
which are then condensed in the extraordinarily compact equation 

V^A = J. (1) 

Although this formulation is valid for Minkowski spacetime, with a different signature 
from 4D0, I shall follow closely that group's exposition. 

Before embarking into the transposition of electromagnetism to 4D0 the paper 
makes a brief introduction to geometric algebra and makes a quick revision of 4D0's 
principles, using the opportunity to express them in the former's formalism. 

2. Introduction to geometric algebra 

The geometric algebra of Euclidean 4-space is generated by the frame of orthonormal 
vectors {cr^}, /i = . . . 3, verifying the relation 

(^ii-CTy = -{cr^ay + a^a^) = 5^y. (2) 
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The algebra is 16-dimensional and is spanned by the basis 

1 scalar 4 vectors 6 bivectors 4 trivectors 1 fourvector 

where / = (Toa"i(72cr3 is also called the pseudoscalar. The elements of this basis are such 
that all vectors and the pseudoscalar square to unity 

{a,r = 1, 7^ = 1; (4) 

and all bivectors and trivectors square to — 1 

Ka^f = -1, = -1. (5) 

It will be convenient to shorten the product of basis vectors with a multi-index compact 
notation; for instance cr^o"iy = cr^,y. 

The geometric product of any two vectors a = a^a^ and b = Wa^ can be decomposed 
into a symmetric part, a scalar called the inner product, and an anti-symmetric part, a 
bivector called the exterior product. 

ab = a-b + aAb, ba = a-b — af\b. (6) 

Reversing the definition one can write internal and exterior products as 

a-b — -{ab -\-ba), aAb — -{ab — ba). (7) 

The exponential of bivectors is especially important and deserves an explanation 
here. If u is a bivector or trivector such that = —1 and 6' is a scalar 

= l-^ + ^----{=cos^^} (8) 

+ u9 — u— + ...{— usm9} 
— cos 9 + usm9. 

Although less important in the present work, the exponential of vectors and fourvector 
can also be defined; if /i is a vector or fourvector such that h'^ — 1 

e- = l + /..4 + /.^4 + ... 

= l + ^ + ^ + ---{=cosh^} (9) 

+ h9 + h— + ...{^h sinh 9} 
= cosh^ -I- hsmh.9. 

The exponential of bivectors is useful for defining rotations; a rotation of vector a by 
angle 9 on the (T12 plane is performed by 
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the tilde denotes reversion, reversing the order of all products. As a check we make 
a = (Ti 



e 



-12e/2^^gai2e/2 _ ^^Qg ^ _ gij^ ^^Qg ^ ^ 



= COS 6' 0"! + sin 6' (72. (11) 

Similarly, if we had made a = (72, the result would have been — sin^ai + cos9a2- 

If we use B to represent a bivector and define its norm by \B\ = (BBY^'^, a general 
rotation in 4-space is represented by the rotor 

R = = cos(|5|/2) - -|- sin(|E|/2). (12) 

\B\ 

The rotation angle is \B\ and it is performed on the plane defined hy B. A rotor is 
defined as a unitary even multivector (a multivector with even grade components only); 
we are particularly interested in rotors with scalar and bivector components. It is more 
general to define a rotation by a plane (bivector) then by an axis (vector) because the 
latter only works in 3D while the former is applicable in any dimension. 

In a general situation the frame may not be orthonormed and we will generally 
define this frame by expressing its vectors in the fiducial frame o"^ 0] 

= h^^aa, (13) 

where /i"^ is called the fiducial tensor. The metric tensor is then defined by the inner 
products of the frame vectors 

g,iu = gt,-gu- (14) 

Complementary we can then define the reciprocal frame by the relation 

g'-g. = 5^. (15) 

This method for frame and metric definition is absolutely general and it is signature 
preserving. Since our fiducial frame has signature (+ + ++) all the frames defined with 
reference to it will preserve that signature. Hestenes 0] discusses curvature in geometric 
algebra formulation but he uses a Minkowski fiducial frame with signature (H ). 



3. Moving frames in 4DO 

In this section we will be looking at the frame applicable to a moving observer and 
making a parallel to Lorentz transformations in special relativity. This is a reformulation 
of the presentation in ^ using the formalism of geometric algebra. We must recall that 
4D0 is characterised by time being evaluated as geodesic arc length in 4-dimensional 
space 

cMt^ = g^^dx^dx". (16) 

From this point onwards we avoid all problems of dimensional homogeneity by 
using normalising factors listed in Table ^ for all units, defined with recourse to the 
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Table 1. Normalising factors for non-dimensional units used in the text 



Length 


Time 


Mass 


Charge 


iGh 


iGh 


jhc 










e 



fundamental constants: Planck constant divided by 27r (K), gravitational constant 
(G), speed of light (c) and proton charge (e). This normalisation defines a system 
of non-dimensional units with important consequences, namely: 1) All the fundamental 
constants, h, G, c, e, become unity; 2) a particle's Compton frequency, defined by 
z/ = mc'^/h, becomes equal to the particle's mass; 3) the frequent term GM/{c^r) is 
simplified to M/r. 

In non-dimensional units equation ()16|) above can be obtained from the equivalent 
vector definition 



As = g^Ax^, (17) 

where ds is the displacement vector and dt^ = ds-ds. Consequently the velocity vector 
is defined by 

V = g^x^, (18) 

where "dot" over a variable means time derivative. 

As an introductory example to moving observer frames we will suppose an observer 

with velocity v = cos6'cto + sin^cxi. Fig. (^. The moving observer is obviously 
stationary in his own frame and so his x° axis must be aligned with v; hence the 
corresponding frame vector must be obtained from v by product with a scalar: go = Xv. 
In 4D0 light is characterised by geodesies with dx° = 0, a necessary condition for the 
3-space velocity to have unitary norm. These displacements must be evaluated similarly 
by the fixed and moving observers so their three frame vectors corresponding to 3-space 
must remain unaltered: (?j = (Xj. 

Any coordinate change cannot alter a displacement so that the latter can be equally 
expressed in either coordinate system 

ds = cr^dx^" = g^.dx''. (19) 

Applying to the case under study 

ds = cTodx" + aidx^ = \{cos9ao + sin 6'cri)dx°. (20) 

As last argument special relativity imposes the invariance of [dx^^ = dt^ — ^(dx*)^; 
in 4D0 we have, from Eq. (jl6|) . dt"^ = (da;'')^-|-^(dx*)^. Consequently, for compatibility 
with special relativity, we must require also that J 

dx° = dx°. (21) 

1 In Almeida 1] a different argument is used to justify this invariance. 
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Figure 1. In his own frame observer O is fixed and go must be aligned with the 
velocity; time intervals must be the same in both frames; dx*^ and time intervals must 
be preserved when the coordinates are transformed; the displacement ^oda::'^ can be 
decomposed into aodx^ + uida;^. 



Combining Eqs. and ()21|1 it must be A = sec 6 and 

90 = o'o + tan 6*0"! = cro(l + tan^^croi). (22) 
The metric tensor components can be evaluated by Eq. (jl4|l 

5foo = 1 + tan^ 9, goi = gw = tan 9. (23) 

The evaluation of dt"^ is unaltered because this is still the fixed observer's time 
measurement performed in the moving observer's coordinates; we call this a coordinate 
change, which is different from a metric change. Using Eq. (fT7j) a general displacement 
is evaluated in the two frames as 

ds = o-odx° + aidx^ = (ctq + tan6'a-i)dx° + (Xidx^; (24) 

and the coordinate conversion is immediate 

dx^ = dx^ -tanMa;°. (25) 

The coordinate transformation for a moving observer seems incompatible with 
the corresponding transformation in special relativity, which is described by a Lorentz 
transformation and produces time dilation. This incompatibility is only apparent since 
the moving observer has no reason to choose a skew frame and evaluates displacements 
with a standard orthonormed frame 

ds = crodx° + cri(dx^ — tan^^dx°) = ds — aitan9dx^ . (26) 
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Notice that the special relativity invariant is preserved in the moving observer's 
evaluation of time 

(ds)!^ = (dt)2 - {dx^f = (dt)2 - (da;^)' = (dx°)^ (27) 

The previous example restrained the analysis to movement along the ui direction, 
described by a velocity vector with ctq and ai components. In general a coordinate 
change defines a skew frame whose vectors result from the transformation applied to 
the element cxo 

gQ = aQ[l + — ] = (Jo = • (28) 

where the bold v represents the 3 spatial velocity components or the velocity vector in 
the non-relativistic sense. The coordinate conversion preserves da;° and for the remaining 
coordinates we have 

k 

dx'^ = dx^-^dx°. (29) 



V 

The interval evaluated by the moving observer in his orthonormed frame is 

V 

-a 

In conclusion, the 4D0 counterpart of a Lorentz transformation in special relativity 
is a two-step process. 



ds = dx^a^, = ds-—. (30) 



• A transformation of the spatial coordinates x*, preserving the fixed observer's 
time measurement; the corresponding frame is called the skew frame and the 
transformation is called a coordinate change. 

• A jump into the orthonormed moving frame with the consequent time dilation. 
This is mathematically a metric change because the length of displacements is not 
preserved. 

A concrete situation where a moving frame had to be considered appeared in 
Almeida 0, Eq. (14)], an equation with the same form of Eq. (j^ . equivalent to a 
frame vector 

9o = cro + (31) 

In this case the movement has a geometric cause and results from the natural (geometric) 
expansion of the Universe when the hypersphere model is assumed. 

4. Electrodynamics as space curvature 

Relativistic dynamics, for cases of isotropic media, is modeled in 4D0 by the frame 

go = noCTQ, gj = Uraj, (32) 

where no and scalar functions of the coordinates called refractive indices. 

Application of Eq. ()16p with c = 1 leads to the time interval 

dt^ = (nodx°)2 + (n,)2 ^(dx^y. (33) 
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If the refractive indices are not functions of x° the geodesies of the space so defined can 
be mapped to the geodesies of the relativistic space defined by the metric 

as is fully demonstrated in Almeida 3:, Sec. 4]. The equivalence between 4D0 and 
GR spaces stops when the metric is non-static, meaning that the refractive indices are 
functions of in 4D0 or t in GR. The two spaces are not equivalent either, for any 
displacements which involve parallel transport. 

The first question we would like to answer in this section is whether a frame of the 
type defined by Eqs. (jH^ is adequate to model the dynamics of a charged particle under 
an electric field; for this we will use the refractive indices 

no = 1 + — , Ur = 1, (35) 
m 

where q is the charge of the moving particle, m is its mass and V is the electric potential, 
including the fine structure constant a. So, for instance, the electric field of a stationary 
particle with charge Q is 

V = ^. (36) 

r 

Notice that the refractive indices are defined for an interaction between two charges, 
the same happening with the space metric. The space that we are defining exists only 
for the interaction under study and it is not a pre-existing arena where the dynamics is 
played. 

The frame of Eqs. ()35j) produces the time definition 

dt' =(^l + ^J {dx^ + 5^(dx^y . (37) 

Following the procedure in Almeida j3| we will find the geodesic equations by first 
dividing both members by dt^ and defining a constant Lagrangian L = 1/2 

1 = 2L = (^1 + {x'Y + ^(x^T. (38) 

Noting that the Lagrangian is independent from there must be a conserved quantity 

l + ^Vx° = l. (39) 

The remaining Euler-Lagrange equations for the geodesies are 

^J = ±(l + 9X.) d,V{xy; (40) 
m \ m J 

Replacing with xp from Eq. 

F = ^(l + '^\ d^V. (41) 
m7"^ y m J 

In the limit of speeds much smaller than the speed of light 7 — 1. As long as qV/m -C 1 
the equation represents the classical dynamics of a particle with mass m and charge q 
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under the electric potential V. Electric potentials typically decrease with 1/r and thus 
the parenthesis can be taken as unity for large distances from electric field sources; for the 
interaction between two electrons, considering the non-dimensional units' normalising 
factors, this condition means distances considerably larger than 2.8 x 10~^^ m. 

Having shown that dynamics under an electric field can be modeled by a suitably 
chosen frame, we need to investigate if the same applies when a magnetic field is present. 
Magnetic fields are originated by moving charges and we have seen how a moving frame 
can be obtained from a stationary one by a transformation applied to its zeroth vector. 
Since electric field dynamics only implies the consideration of refractive index no, the 
other refractive index Ur remaining unity, it is natural to admit that an electromagnetic 
interaction could be modeled by the following frame 

5-0 = o-Q H , = CTj, (42) 

m 

where A = A'^a^ represents the vector potential; note that the vector potential is referred 
to the orthonormed frame cr^ and not to the skew frame Qf^. The equation above is 
equivalent to the definition of a fiducial tensor but the vector potential approach is 
more convenient for the derivations that follow. 

The metric tensor elements are obtained, as usual, by Eq. (jHj) 

/ qA^Y qA^ ^ , , 

5'oo = IH , 9^3 = 9jo = , 9jk = Ojk- (43) 

\ m J m 

Instead of using the geodesic Lagrangian to find its Euler-Lagrange equations, we will 
follow a different procedure, starting with the velocity vector from Eq. p8|) with the 
frame vectors from Eqs. ()42|1 

= ^ i;0 + a.x''. (44) 
m 

We will now define the vector derivative 

V = a^d^. (45) 
This is used for the derivation of an identity relative to acceleration 

V = v^a^ 

= ^,v^'v^[a,■{a,^a^)] (46) 

= v-iyAv). 

Applying to Eq. 



VA( + 

m 



(47) 



The exterior product with the second term inside the parenthesis is necessarily null 
because it implies deriving the coordinates with respect to other coordinates and not to 
themselves; so we have finally 

iJ='^v.{y^A). (48) 
m 
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It is now convenient to define the Faraday bivector F = V /\A. F is necessarily a 
bivector because it is defined as an exterior product of two vectors; we can separate it 
into electric and magnetic components 

F = E^ajo + BHak. (49) 

E^ and are the components of electric and magnetic field vectors, respectively; 
i = c"i23 is a special trivector incorporating the 3 spatial frame vectors, so that iak is 
a spatial bivector where the component cr^ is not present. Using boldface letters to 
represent 3-vectors, the Faraday bivector can be rewritten as 

F = Edo + iB. (50) 

And the acceleration becomes 

ij = ^y.(Eao + iB). (51) 
m 

This form of the Lorentz acceleration is relativistic because the norm of the velocity 
vector will be kept equal to unity at all times. 



5. Maxwell's equations 

Using the vector derivative defined in Eq. ()45|) we write the derivative of Faraday bivector 



VF = V-F + VAF, (52) 

The internal product in the second member is expanded as 

V-F = dkE'ao - doE'ak + {d^B^ - d2B'')a, (53) 

+ {diB^ - d-iB^)a2 + {d2B^ - diB^)a^. (54) 



A careful look at the second member shows that the first term is the divergence of E 
multiplied by o"o. The second term is — SqE, but we must take into account that in 
a stationary frame dt = dx°, so this term can be taken as the negative of the time 
derivative in such frame. The last 3 terms represent the 3-dimensional cross product, 
for which we will use the symbol "x". If the bold symbol "V" represents the usual 
3-dimensional nabla operator, the equation can be written 

V-F = V-E(To + Vx B-(9oE. (55) 

Defining the current vector J = paQ + J, for a stationary frame, the first two Maxwell's 
equations are expressed by 

V-F = J. (56) 
For the remaining two equations we have to look at the exterior product VAF 

VAF = {d,E^ - d2E^)ao23 + {diE'' - d,E')aou (57) 
+ (^2^1 - d,E'')aou + idkB'' + doB'^Iak (58) 
= i(VxEao + V-B + (9oB(To). (59) 
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The remaining two equations are expressed as VAF = and the four equations can 
then be condensed in the very compact form 

VF = J. (60) 

Additionally, it is possible to choose A such that V-A = 0, so that an equivalent form 
of the four equations is 

V^A = J, (V-A = 0.) (61) 

Solutions of special interest are those where the second member is null (J = 0). 
Some solutions can be found easily by splitting V^A in its 0th and spatial components 

W^A = dooA + V^A, (62) 

where the second term is a standard 3D Laplacian. One possible solution arises 
immediately when OoqA = uJ^A with uj real; 

V^A = -u?A, (63) 

results in a well known Helmholtz equation. Looking for a particular solution we assume 
only and x^ dependence to find 

A = Aoe-"^^"-""^'. (64) 

Where = —\ and u must commute with Oj,. The two possibilities are either u = i or 
u = ia^; we will choose the latter. Since we know that A is a vector it must be 

Aq = aai + (3a2-, (65) 

an eventual component aligned with ctq has been ignored. 

In order to make interpretation of Eq. easier we shall now write it in a slightly 
modified form 

A = e-'""°e'"^""'/2(a(Ti + /5a2)e-'"^""'/l (66) 

It is now apparent that the first exponential factor is evanescent in the positive 
direction and grows to infinity in the negative direction; this is an uncomfortable 
situation that will be resolved below. The remaining factors represent a vector that 
rotates in the ia^ plane, along the x^ direction, with angular frequency uj. With an 
adequate choice for the origin of x^ it is possible to make (3 = 0. 

Before we proceed to the analysis of evanescence along x° let us consider consider 
the expression for A when x° is constant; setting a' = a exp{—u!x^) it is 

A = aVie"'"^""'. (67) 

Along the direction normal to the rotation plane we can define time to be equal to the 
distance traveled, t = x^ and so A can also be expressed as a time function 

A = a' aie-'"^'^'. (68) 

In spacetime formulation time is a coordinate and we express the fact that t = x^ hy 
giving A opposite dependencies in the two variables 

A = aVie-^"3'"(^'-*). (69) 
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This is the spacetime formulation for a circularly polarised electromagnetic wave j^. 

In Almeida the Universe is modeled as an hypersphere whose radius is the 
coordinate . For displacements that are small compared to the hypersphere radius 
it is possible to use length coordinates rather than angles, resulting in the frame 
vector of Eq. (p?T|) . A comparison between this equation and Eq. ()42|1 shows that the 
Universe's expansion is the source of a vector potential A = x'^ak/x^; the derivations 
in the previous section show that there are solutions to Maxwell's equations that force 
the vector potential to rotate and the frame vector go with it. These solutions force 
x° to stay constant and so they are constrained to hyperspherical surfaces of constant 
radius. This solves the unanswered question raised in Almeida whereby photons were 
constrained to great circles on the hypersphere, rather than to geodesic straight lines. 

We can now turn our attention to the evanescence problem or rather to the problem 
of an amplitude growing exponentially in the negative x° direction. This is the direction 
pointing towards the hypersphere centre and suggests that the approximations made 
in order to use flat space are the source of the inconsistency. Recalling that in the 
hypersphere model the Universe is expanding at the speed of light, we can expect that 
there are resonating modes inside the hypersphere, which are evanescent to its outside 
due to continuity on the hypersphere border. A correct solution of Maxwell's equations 
should not ignore this fact and will be pursued in forthcoming work. The resonating 
modes are so closely spaced with the current size of the Universe that we can actually 
consider a continuous distribution, but they manifest themselves as electromagnetic 
waves when space is artificially flattened. 

6. Conclusion 

The formalism of geometric algebra was used in this paper as useful mathematical 
tool for writing complex equations in a compact form and facilitate their geometrical 
interpretation. Although most readers will not be familiar with this algebra, we think 
that the simplifications achieved through its use are well worth the extra effort of 
learning some of its rudiments. A reformulation of Maxwell's equations in 4-dimensional 
Euclidean space is the goal of this work but an introduction to the discipline of 4- 
dimensional optics was felt necessary, namely to explain the new transformation for 
moving frames and its relation to Lorentz transformations in special relativity. 

Through the consideration of an electromagnetic vector potential embedded in one 
of the frame vectors, the paper shows that relativistic electrodynamics can be derived 
from curved space geodesies. In this way the Lorentz force acquires the characteristic of 
an inertial force, which is best described by inertial movement in curved space. Maxwell's 
equations are then established in this space, with geometric algebra achieving the feat of 
their condensation in the equation V'^A = J. Solving the equations for free space leads 
to rotation of frame vector go on a plane lying on 3-space, with rotation progressing at 
arbitrary frequency along the direction normal to the rotation plane. These solutions 
are evanescent on the positive x^ direction and grow to infinity in the opposite direction. 
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The inconsistency is attributed to artificial flattening of the hypersphere space proposed 
in previous work |^ and further work is suggested to fully clarify this point. 
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